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We show that describing the screened electrostatic interactions in a periodic stack of rigid parallel 
identical charged plate-like colloids within a local density functional theory approach, generically 
leads to a swelling behaviour. Within the same framework, we find that the effective potential 
between a pair of such plates immersed in an electrolyte is repulsive at all distances. This general 
result is in contradiction with a theory put forward by Sogami, Shinohara and Smalley, that we 
criticize and thoroughly reconsider to show how the two approaches can be reconciled. 



I. INTRODUCTION 



The counter-intuitive phenomenon of effective attractions between like-charges immersed in an electrolyte has 
recently attracted considerable attention, and questions one of the fundamental tenets of colloid science jl] g. The 
understanding of the underlying mechanism is crucial for a correct description of a vast variety of industrial and natural 
processes, in particular the stability and phase behaviour of colloidal dispersions. In this article, we investigate the 
stability of electrostatically swollen stacks of regularly spaced planar colloids or membranes of infinite lateral extension 
in an electrolyte. Such a one-dimensional geometry describes the early stages of clay swelling |j| but is also relevant for 
lamellar phases of charged bilayers Q . At the simplest level of mean- field description, with neglect of excluded volume 
and Coulomb correlation effects, the inhomogeneous density profiles of microscopic co- and counter-ions forming the 
electric double layers around the planar macroions can be obtained analytically |8[ . The resulting Poisson-Boltzmann 
osmotic (or disjoining) pressure in the stack, is invariably found positive qJ, which is the signature of the tendency 
to swell. This behaviour is reminiscent of the repulsive nature of pair interactions (as we shall see below, the two 
phenomena are however distinct), established within the same mean- field framework and recently extended 

beyond mean-field Jl2| , with the possible inclusion of approximate correlation contributions. Note that more refined 
incorporation of discrete solvent effects by adding a bilinear non local term to the free energy of standard Poisson- 
Boltzmann theory has shown the possibility of a net attraction (negative pressure) at high surface charges Jl3| . 

The purpose of the present investigation is twofold. First, we apply the general local density functional formalism 
of Ref. (l2| to provide a prediction for the stability of stacked or lamellar phases. We show in section || that the 
interactions of electric double layers generically lead to swelling whereas describing an equilibrium spacing between the 
membranes either requires the inclusion of non electrostatic forces (usually van-der-Waals like) or the consideration of 
more refined theories. Second, the pair potential problem is considered in section III with a resulting effective repulsion 
at all distances, within the same framework that encompasses in particular the non-linear Poisson-Boltzmann theory. 
This statement is consequently in contradiction with the attractive behaviour reported by Sogami, Shinohara and 
Smalley ■ We show that their controversial finding is attributable to a confusion in the thermodynamic potential 

describing the electrostatic situation under study. 



II. SWELLING OF A LAMELLAR STACK 



Before setting the framework of the analysis, it is worthwhile to point out that the results obtained in 10|-|l^] for a 
pair of colloids do not apply to the regular stacked situation under scrutiny here (the precise definition of the geometry 
is given at the beginning of section IIA). First of all, the mechanical route followed in Refs. fic|-|l^| (deriving the 
interactions by integration of the stress tensor over the colloids' surfaces) would yield a vanishing effective force in a 
regular stack, with cancelling contributions acting on both sides of the platelets, and is thus uninformative about the 
stability of the array. Moreover, the situation considered in [l^] is that of two colloids immersed in an electrolyte, 
confined in a cylinder with an axis parallel to the colloids' line of centers, and of infinite extension along this axis 
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[see Figure Q-a)]. The assumption of infinite length is crucial for the validity of the effective repulsion (see Appendix 
|A|) and the corresponding confined two-body problem does not include the multi-body stack geometry, for which 
the confining cylinder (Wigner-Seitz cell) would be a slab of finite length. This point is in contradiction with one 
erroneous conclusion reached in Ref . [Elf . 



A. Density functional theory formalism 

We consider a regular succession of rigid infinite parallel plates located at z = 2nh (n G N). Each plate carries a 
uniform surface charge a, and releases its counter-ions in the electrolyte solution, considered to be a mixture of N 
microions, where species a has charge number e a and local density n a (z). We write the free energy p6[ of the total 
charge distribution in the Wigner-Seitz cell around a given platelet (e.g. that situated at z — for which the cell is a 
slab extending between —h and h) as 



1 



f[{n a (z)}} dz + - / p c {z)^{z)dz 



(1) 



-h 



where p c (z) is the local total charge distribution, including the microions and the plate (with global electroneutrality 
j, p c = 0). J- is expressed per unit area of the macroscopic surfaces. Within the primitive model of electrolytes p2| , 
whereby the solvent is regarded as a mere dielectric continuum of permittivity e, the total electrostatic potential if(z) 
is the solution of Poisson's equation 



V 2 V> = Pc{z) = 

£ £ 



aS(z) + e a n a (z) 



(2) 



and can be cast in the form ij){z) = j , p c (z') G(z, z') dz' where G is the appropriate Green's function. The boundary 
condition of vanishing electric field at z = ±/i is fulfilled by d z ip. This mean-field like reduction of the problem 
materialized by the introduction of the Wigner-Seitz slab, may be corrected by correlation terms included in the free 
energy density /. Moreover, even if the van der Waals energy term J p c ip = J p c Gp c is of mean-field form, correlation 
terms can be accounted for provided they translate into a local correction to the free energy, in the spirit of the 
approaches depicted in Refs ||| . Consequently, the term J / in Equation (|l|) does not reduce in general to the 
entropic microions' contribution (as in Poisson-Boltzmann theory |^| ) but may also include both non mean-field 
energetic and entropic effects. Of course, the present formalism encompasses the non-linear Poisson-Boltzmann and 
Modified Poisson-Boltzmann theories |B7LE8fl. 



B. Osmotic pressure 



For both canonical (fixed ionic content of the electrolyte solution) and semi grand-canonical (when the solution is in 
osmotic equilibrium with a salt reservoir) descriptions, the optimal density profiles n* a {z) are obtained by minimizing 
the functional 

n({n a }) = T({n a }) - J2^ N ^ \N a = J ^n a dz^j, (3) 

where p a is either the Lagrange multiplier associated with the constraint of fixed average concentration, or the 
chemical potential of species a in the reservoir. The resulting stationary equations read for all species a 

df 

— 1- e a '4>(z) = p a . (4) 

These implicit relations between the electrostatic potential and the densities {n a } allow to solve Poisson's equation 
(H) and compute the optimal profile n* a (z) (hereafter noted n a (z) without ambiguity), from which we deduce the 
Helmholtz free energy F(h, a, T, {^V Q }) = T({n* a }), with F = U — TS {U is the total internal energy, T the absolute 
temperature and S the entropy of the total micro- ions charge distribution in the cell). 

From the knowledge of F, the definition of the thermodynamic potential 1Z (which variations correspond to the 
reversible work performed by an operator and thus define the osmotic pressure) requires the specification of the 
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thermodynamic situation under consideration. In the canonical case with constant charge plates, TZ — F p9[ . If 
on the other hand, the platelets are held at constant potential (e.g. by an external generator), TZ is the Legendre 
transform of F with respect to the surface charge P9,pQ| , namely 



OF 

TZ = F(a,h,T)~a— 
do 



F(a,h,T) — ai\) v (cr, h,T), 



(5) 



h.T 



where ip v = tp(z = 0) denotes the surface potential. In the opposite case of a system in equilibrium with a salt 
reservoir, 



TZ = TZ a = fl({n a }) = F — fi a / n a (z) dz 
TZ.,p = F - J^Ma / n a (z)dz 

J-h 



TZ = 



In any case, the osmotic pressure is defined as 



dF 



at constant charge a 
at constant potential 



n 



2 ~dh' 



(6) 
(7) 

(8) 



where TZ depends on the electrostatic situation considered, as explained above. It is however important to stress that 
the osmotic pressure should not depend on the electrostatic situation under study, as becomes clear below. Theories 
that do not result in the above fundamental invariance of II can consequently be disposed of. 

The free energy variation induced by changing the inter-membrane distance is computed in Appendix |b] with the 
result 



SF = 

so that the osmotic pressure reads 



2Sh + ^2fi a SN a + ip p 5a, 



(9) 



- z—h 



n 



df_ 



(10) 



independently of the situation of constant potential or constant charge considered. Not surprisingly, the pressure is 
related to the Legendre transform of the free energy density, as for ordinary homogeneous gases |30| and equals the 
local stress ir(z) at the mid-plane between the membranes: 



d f 

II =ir(z = h) where n(z) = -f[{n a }(z)] + 2J n a (z)— —(z) 



dn a 



(11) 



C. Stability analysis 



When the stack is in osmotic equilibrium with a salt reservoir, the comparison of II obtained in section II B with 
the pressure exerted by the reservoir quantifies the tendency towards swelling (II > II rcs ) or collapse (IT < II res ). 
For consistency, the (neutral) reservoir needs to be described within the same framework as the electrolyte around 
the platelet. The remainder of this section is devoted to the proof that LT is extremal in the reservoir, and that this 
extremum is a minimum under fairly general conditions. We first analyze the ip dependence of the pressure tt, defined 
by n(ip) = tt(z) where ?/; = ^(z) is the solution of Poisson's equation (g). From Eq. ( [Tl| ) we have 

that can be recast making use of the stationary condition (Q) 
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()tt 



N 



-^n a e a (13) 

a=l 

-Pc(z) outside the plate (z ^ 0). (14) 



7r thus goes through an extremum in the reservoir (II rcs by definition) where the charge density vanishes (unlike at 
the mid-plane z = h where p c ^= 0). Relation (p~4]) together with Poisson's equation implies that the local stress 
introduced in ([ll]) fulfills the mechanical equilibrium condition for the fluid of microions: 

*(*)-£(^Y = n, (is) 



8tt \dz 

which is the simple form taken in a one-dimensional problem by the condition of vanishing divergence for the gener- 
alized Maxwell stress tensor in the dielectric medium 



The second derivative of ir can be obtained by introducing the Hessian matrix N x N 



that obeys the relation 



d 2 f 

H a , = —-j— (16) 



/3=1 y 



P-^T = ~ ea ( 17 ) 



as can be seen from Eq. m). We thus have 



d n _ ST h dria dn P n»l 

which is a positive quantity when the matrix of second derivatives H a p is positive definite. Under this assumption, 7r 
is a convex-up function of the potential and the extremum attained in the reservoir is therefore a minimum, so that 

n > n ros . 

In spite of the differences outlined at the beginning of section |n] between the pair potential and the multi-body 
stacked problems, the method employed here is close to that of Ref . , and we shall also distinguish two different 
situations to conclude with the stability analysis 

• Case ah the free energy density does not depend on the elementary charge e (as in mere mean- field treatments 
P|, p2|j27|j2^ | ) . In the limit e —> (at fixed valency e Q /e), we obtain a locally neutral mixture where the total 
free energy functional reduces to the contribution J f. The thermodynamic stability criterion of this neutral 
mixture implies the positive definitencss of the matrix H a p that is independent on e. 

• Case b): correlation or fluctuation effects are taken into account with a resulting e-dependent free energy density. 
It is no longer possible to find an uncharged mixture of micro-species described by the same density /. The 
thermodynamic stability condition of the full functional (Q) involves the Green's function G(z,z') and does not 
imply the positive definiteness of the Hessian H. However, the convexity of / with respect to density variations 
is generally fulfilled by the approaches proposed in the literature, either in the full density range |25| or for the 
small plasma coupling parameters relevant to colloidal dispersions 

From the above discussion and the positive definiteness of the Hessian matrix H , we conclude that the stack generically 
exhibits a tendency towards swelling. 



III. ADIABATIC PAIR POTENTIAL 



In this section, we consider within the density functional formalism of Eq. (|l|), the problem of the effective 
interactions between a pair of parallel charged plates immersed in an electrolyte solution of infinite volume (no 
confinement). The two rigid plates with distance 2h are supposed to be of negligible thickness and uniform surface 
charge (as in section Hi); they divide the electrolyte solution into two disconnected regions (inner region with \z\ < h 
and outer region with]z| > h), which does not correspond to the situation analyzed in Refs. |10|-p^|. However, the 
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generalization to the present case is straightforward and yields an effective repulsion as long as both regions are in 
contact with salt reservoirs imposing the sa me ch emical potential for micro-species. Indeed, from the computation of 



the free energy variation detailed in section II B , the effective force can be written 



From the vanishing of the electric field at z — and \z\ — ► oo, we have IIi n = n(z — 0) and n out = it(\z\ — > oo). For 
\z\ — » oo, the charge density vanishes so that n out equals the osmotic pressure in the salt reservoir under consideration. 



The argument of section [I C indicates that IIi n > n out and that the interactions are repulsive under the assumption 
of positive definiteness for the stability matrix d^of. This last condition is obeyed by Poisson-Boltzmann (PB) theory 
independently of the valency of the microions [see Eq. ( |20| ) below]. Our result is consequently in contradiction 
with the "long-range weak attractive part of the free energy" reported in Jl4|,|l5| for the same system treated at the 
level of Poisson-Boltzmann. The work of Sogami et al. has already been criticized p6|^I^], but in our opinion, the 
subsequent controversy pfj| dwells on ambiguities on the thermodynamic potential that should be considered, which 
to our knowledge, have not been explicitly pointed out so far. It thus seems worthwhile to restrict to PB theory 



and devote the remainder of this article to briefly revisit the model of [|14 15], introduced to describe the swelling 
behaviour of n-butyl- ammonium vermiculite gels. 

Within PB mean-field theory, the microions are considered as an ideal gas and density fluctuations discarded, so 
that the free energy density does not include any correlation term and reduces to the entropy of an ideal mixture: 

JV 

f({n a }) = p- 1 E n « t ln K A ') - !] ' ( 2 °) 

a=l 

where the (irrelevant) lengths {A Q } involve the masses of microions and (3 — l/(kT) is the inverse temperature. The 
stationary condition (|]) translates into 

n a (z) = n° a exp(-/3e Q -0) (21) 

and the local osmotic stress in ( |ll| ) is given by the ideal equation of state ir(z) = kT^2 a n a (z). The electrostatic 
potential ip is chosen to vanish for \z\ — ► oo, so that n° is the density of species a far from the plates with a 
corresponding chemical potential 

fi a =kTln(n° a Ai). (22) 

As a result of the global electroneutrality condition, it can be checked that the final free energy [expression ( p3|) below] 
is independent on an arbitrary shift of the potential ip [in which case the densities appearing in Eq. ( |22[ ) are simply 
the prefactors of the exponentials in (^lj)]. It is convenient to use the relation between the local ionic densities and 
the electrostatic potential to recast the Hclmholtz free energy F = U — TS in the form: 

F = HM) = ai> v - J (j^J dz + Y,N a (p a - kT) . 

/OO 
Ildz + ^2N a {fi a -2kT), (23) 
-oo a 

where Eq. ( |i~5| ) and the ideal equation of state for 7r(z) have been used in going from the first to second line. Note 
that it is understood that the osmotic term II in ( p3[ ) takes the value IIi n (resp. n out ) for \z\ < h (resp. \z\ > h). 
Strictly speaking, expression (^) diverges (if salt is added to the electrolyte, some of the quantities N a are extensive 
with system size) . This feature can be circumvented by computing the excess free energy with respect to a well chosen 
reference system (for instance the system with same bulk densities in the absence of the plates). Once the solution 
of Poisson's equation is known (see || for a review of the standard solutions, including the present geometry), F 
may be computed from equation (^3|). Following this route, we readily recover the Helmholtz free energy obtained by 
Sogami et al. by means of a charging process. In Equation (p3|), a = a m + <r out stands for the total surface charge on 
a platelet, including both the different inner (cri n from z — h~) and outer (<7 ut from z = h + ) contributions, denoted 
respectively Zi and Z Q in Jl^JTo] ]. 

In the model of |M|]l~5| ], the electrostatic potential is imposed to be continuous throughout the system and the 
surface potential %p v is independent on the distance 2h between the plates. Moreover, the numbers of microions 
between the plates also depend on h with a fixed chemical potential, given by (K3) where the n° a can be considered as 
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the ionic densities of a reservoir in chemical equilibrium with both the inner and outerparts of the electrolyte solution 



erps 

around the plates. Consequently, the thermodynamic potential TZ^ defined in Eq. (|7j) should be used in computing 
the force whereas F has been considered in |l4|,[l5| . The definition of the model imposes that when h changes at fixed 
ip v , Z a is constant while Zi varies, so that the two situations of constant charge and constant potential are equivalent 
for the outer part of the system (\z\ > h), but not for the inner part. The reversible work performed by an operator 
changing the distance 2h between the charged plates is therefore given by the variations of 
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■4, 



F\r\ ~\~ Fry 



E 



Ma 



(24) 



whereas F — F in + F out has been considered in [ [1411151 . Omission of the chemical potential terms fi a n a in the right 
hand-side of Eq. ( p4| ) leads the force computed in ||14| , ^5| to depend on the masses (through the lengths A Q ), which is 
impossible in an equilibrium statistical mechanics theory, as already pointed out in |H| . With the aid of the relations 
given m @ [e.g. Eqs (63), (64) and (65) valid in the case of counter- ion dominance between the plates] , an explicit 
computation of the force F z = —d1Z^,/d(2h) yields the standard expression 



F z = ir(h~) — — (-jr-) 



* + >-s(£L. 



8^ 



dy - 

dz 



=h+ 



dz . 



-h- 



(25) 
(26) 



where the continuity of the co- and counter-ion charge density across the membrane (resulting from the imposed 
continuity of the potential) has been used in going from (|25|) to (|26|). The compatibility of Eq. (|25| ) with ([11]) is 
transparent. In terms of the variables used in [ p^l5| ], where $(0) denotes the reduced potential at mid-distance 
between the plates, we get 



F z = -eHzl 



Zf) = 4n fcTsinh 2 



$(0) 



>o, 



(27) 



hence an effective repulsion at all distances, that is in the present case entirely due to the electrostatic pressure (the 
osmotic contribution cancelling on both sides of the plates). 

It can be checked that in the dual situation where both the inner a- m are outer <r ou t surface charges are held constant, 
the effective force is still given by ([25|). This translates into the Legendre identity 



dh 



!<7out 



d_ 

dh 



- ip. 



'<P V : ff OUt. 



8 (7^/, + (TinVV ) 



dh 



> CT out 



(28) 



^in > CT out 



that can be considered as a test for the consistency of the thermodynamic potential used. In the limit of counter-ion 
dominance between the plates investigated in (TJ] , the identity ((2^) can be checked explicitly with the thermodynamic 
potential IZ^i, used here. 



IV. CONCLUSION 



Describing the interactions between the electric double-layers around charged planar colloids with a local density 
functional theory for ionic screening and the primitive model of electrolytes, we have shown that a regular stack of 
such plates generally displays a swelling behaviour when electrostatic forces alone are taken into account. Within the 
same framework, the effective pair potential is found to be repulsive at all distances, in contradistinction with the 
results derived in |l4|,|l5|]. The repulsive interactions are evidenced without resort to an explicit solution of Poisson's 
equation, and are related to the convexity of the underlying free energy functional, as already noted in JI4] , 

Of course, the present result does not preclude the possibility of effective attractive pair potentials. In particular, 
a drawback of the theories encompassed by Eq. ([!]) is that the direct correlation function c^(ri, r^) defined as the 
second functional derivative of T is necessarily mean-field like: 



.(2) 

-a,/3 



(ri,r 2 ) = -- 



d 2 f 



-<5(ri-r 2 ) - e Q e/3G(ri,r 2 ), 



(29) 



dn a dnp 

with a ^-correlated short range part. This shortcoming may be circumvented by the inclusion of non-local terms in 
the theory, e.g. in the spirit of the weighted density approximation J33|,|34j]. However, we expect the local formalism 
considered here to be instructive for inter-plate separations much larger than the ionic size. 
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APPENDIX A: 



References lid proved that the effective interactions between a pair of like-charged colloids immersed with 
counter-ions and salt in a confining cylinder of infinite extension [see Figure [lj-a)] were repulsive. We show here that 
the result does not hold for a finite length cylinder. 

We consider two colloids confined in the cylinder of length 2h represented in Figure |l|-b), in the specific case of 
Neumann boundary conditions on the surface £ (vanishing normal electric field). Due to the mirror symmetry of the 
charge distribution, the electric field has no z-component in the plane z = 0, and the problem is equivalent to that of 
a unique colloid (say colloid 1) in a sub-cell cylinder of length h, with again Neumann boundary conditions [the right 
half part of the cylinder, indicated by a dotted rectangle in Fig. |l|-b)]. We assume the effective force acting on colloid 
f repulsive (Ff > 0). Then, wc translate along the z axis the left half cell {—h < z < 0) by a distance 2h, keeping 
the right half fixed. The electrostatic potential around colloid 1 is unaffected, so that the effective force is unchanged, 
corresponding now to an effective attraction. Any repulsive configuration with Neumann boundary conditions on the 
confining cylinder can then be mapped onto an attractive one (from the construction of the mapping, it appears that 
this feature disappears in the limit h — » oo). The hypothesis of infinite length is thus a key ingredient of the proofs 
in Refs. ]Tc|-|l^|. It is moreover worthwhile to note that numerical solutions of the non-linear Poisson-Boltzmann 
equation for finite-size disc-like clay platelets confined in a finite length Wigner-Seitz cylinder, show repulsive effective 
pair forces when the distance between the clay particles is smaller than the half length h of the cylinder [B5| . 




FIG. 1. Illustrative side view of the cell models considered. Mirror symmetry with respect to the mid-plane z — between 
the two colloids is enforced (this requirement is an important ingredient in the proofs [|I(fl-Pl2|). 



APPENDIX B: 



It will be shown that upon modifying the distance between the plates (h — > h + Sh), the Hclmholtz free energy 
changes according to Eq. (|9|) . The present derivation bears some similarities with other ones in the related context of 
Poisson-Boltzmann cell theory |32|. From the definition of the free energy, Equation (|l|) and the symmetry z <-» —z, 
we have 



h ph+Sh ^ r^ 1 ph-\-5h 



SF = / 5fdz + 2 fdz+- / 5( Pc ip)dz+ / p^dz. (Bf) 

J-h Jh 2 J_ h J h 

Consider first the energetic contribution. Making use of Poisson's equation (||), two integrations by parts [with 
vanishing terms (<9 z V0z=±fJ yield 



h ph ph+Sh 

5(p c ip) dz= / ipSp c dz + / pcipdz (B2) 



h Jh 



/h N ^ „h+Sh / df \ 

ipSpcdz + X] y n ° y 1 " _ dz ' ( B3 - ) 
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The stationary condition (H) was used in going from (B2) to (B3). Similarly, with N a = J_ h n a dz and since p c 
reduces to the microions charge distribution ^2 a e a n a outside the platelets (in particular between h and h + 5h) 



Sfdz 



N „h N / rh+Sh 

^ / (fJ-a- e a ip) 5n a dz = ^ fi a I 5N a - 2 1 

a=l J ~ h a=l \ Jh 



i)5p c dz + ^(z = 0)8cr. (B4) 



Gathering results, we obtain equation 
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